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Abstract 



A neutral dark matter particle may possess an electric dipole moment (EDM) or a mag- 
netic dipole moment (MDM), so that its scattering with nuclei is governed by electromagnetic 
interactions. If the moments are associated with relevant operators of dimension-5, they may 
be detectable in direct search experiments. We calculate complete expressions of the scattering 
cross sections and the recoil energy spectra for dark matter with these attributes. We also pro- 

m 

£T) • (CFF) which is related to the charge radius defined by an effective dimension-6 operator. We 



show that a 7 GeV dark matter particle with an EDM, MDM or CFF easily reproduces the 
f^ ■ CoGeNT excess while remaining consistent with null searches. 



1 Introduction 

The nature of the dark matter (DM) particle is unknown. There are many well studied scenarios 
guided by theoretical niceties such as the Minimal Supersymmetric Standard Model, models with 
extra dimensions, and the Little Higgs Models with T-parity, all of which posit that the DM particle 
interacts primarily via weak interactions. On the experimental side, direct DM searches are looking 
for signals of the recoiling nuclei from DM-nucleus scattering. We take the uncommon view that the 
scattering process may be electromagnetic in nature. The interaction occurs through the electric 
dipole moment (EDM) or the magnetic dipole moment (MDM) of dark matter. The interactions 
are described by dimension-5 operators for non-self-conjugate particles, such as Dirac DM, but not 
for Majorana DM. The EDM and MDM of DM can be induced by underlying short distance physics 
at the one-loop order or higher. As there are no strong reasons against large CP violation in the 
DM sector, we cannot neglect the EDM possibility. In fact, for comparable short distance cutoffs, 
we find that an EDM may give the dominant contribution to DM-nucleus scattering because it 
directly couples to the nuclear charge. In the EDM case, the recoil energy Er distribution is highly 
enhanced in the low Er region as 1/{v 2 Eji) (where v r is the speed of the DM particle in the rest 
frame of the nucleus) in contrast with that from MDM which goes as 1/Er. 

While DM-nucleus scattering due to the MDM of DM has been studied extensively in the 
literature, the relevant formulae with the correct dependence on the nuclear charge Z and nuclear 
moment Hz, A are not available. In this work, we provide analytic expressions for the scattering cross 
sections for the MDM and EDM cases with careful expansions in the relative velocity and recoil 
energy. To complete the treatment of electromagnetic properties of DM, we extend our analysis to 
the dimension-6 operator, which is the electric charge form factor (CFF) slope or the charge radius 
of the neutral DM. The operator has a structure that is similar to that of the spin-independent 
(SI) interaction. 

2 Electric dipole moment of dark matter 

The effective non-relativistic Hamiltonian of the EDM of a particle with spin S is 

H cS = -dE-S/S, 

with the normalization chosen to agree with the standard form for a spin ^ particle, i.e., —dE ■ a, 
where d is dimensionful, and the electromagnetic energy density is (E 2 + B 2 )/2. As the electric field 
of the nucleus E = — grad (f>, we identify the gradient with the momentum transfer q. Therefore, 
the direct scattering of the DM particle % and the nucleus N, %Af — > %Af, via the interaction 



between the DM electric dipole (with moment d x ) and the nuclear charge Ze (with e 2 = 4ira) is 

M = S-iq (jg^Ze 

It is important to note that the momentum transfer q is Galilean invariant and is conveniently 
related to the center-of-mass momentum p. Hence dq 2 = 2\p\ 2 dcos6 = 4\p\ 2 ^, and \p\ = m r v r , 
where m r = ™ Arn ^ is the reduced mass; ijia and m v are the masses of the nucleus and DM 
particle, respectively, q 2 ranges from to (2m r v r ) 2 , as is easily checked in the center-of-mass frame 
in which the momenta have equal magnitudes m r rvLJ We use the spin relation, 

TriPS*) = (2S + l)6 ij S(S + l)/3 , 

to obtain the spin-averaged differential cross section, 

d*E DM ( X N) = ^ 2 e 2 ^^\G E {q 2 )\ 2 . (1) 

We have included the nuclear charge form factor \GE{q 2 )\ 2 to incorporate elastic scattering effects 



off a heavy nucleus 



Accounting for a difference in convention in the definition of e, our formula 



agrees with that of Ref . [6] . 



1 The DM particle velocity in the frame of the galactic halo is usually described U by a Maxwellian distribution, 

/ (v)dv = - 7 =- s e <°dv, 

where the most probable velocity is vo — 230 km/s, and the distribution is cut off at the escape velocity 
Vase — 600 km/s. The normalization TV is close to 1, or more precisely, 

7V= I . 

erf(^ csc /-u ) - ■^(v csc /v )exp(-vZ sc /v$) 

The one- variable velocity distribution is 

G 4v 2 N - v i/ v 2 

fx(v)dv — „ ,_ e °dv. 

"ov^ 

However, since the solar system is moving at a speed ve — 244 km/s with respect to the halo [2] and ve ~ Vo, we 
need to use a more relevant distribution of the relative velocity, v r — v — ve- As an approximation we ignore the 
seasonal motion of the Earth around the Sun with a relative speed of about 30 km/s. Then, 

f(v r )d v r = / (v)d v = / (v r + VE)d v r because d v = d v r . 

For a fixed v r , we integrate the polar angle between v r and ve to obtain 

fl(Vr)dV r = N VrdVr { e -(min(v r -v E ,v csc )) 2 /v% _ e -(min(v r + v E ,v CBC )f /v%\ 

ff v r > v csc + v E , ,fi(v r ) — 0; see Refs. [U[3]. 

2 A good nuclear form factor can be found in Ref. @j. The spatial charge distribution is parameterized by the 

Fermi distribution p(r) = po/(l + e i - r ~ c ^ a °), where the radius at which the density is po/2 is c = (1.1&4.3 — 0.48) fm 

and the edge thickness parameter is a^ = 0.57 fm [5]. The form factor that is valid for nuclei with a well-developed 



We relate q 2 to the nuclear recoil energy in the lab frame, q 2 = 2mAEji 1 and find 
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The \/(y 2 Ept) dependence is characteristic of the EDM of the DM particle 

To have an EDM, the DM particle cannot be self-conjugate. Consequently, for S 
particle has to be Dirac. Note that the spin factor ^±7- becomes 1 for 5 = \ in our numerical 
illustrations. Our result also applies to the anti-dark matter particle under the assumption that 
CPT is conserved. 



core (i.e., with atomic masses above 20) is obtained by the Fourier transform in the limit c ~S> ao, 



G E (q) = 



nao sin(qc) cosh(7raoq) cos(qc) 
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Note that G E (0) = 1. 

3 The differential reaction rate (per unit detector mass) is 



dR po 1 f°° . , > da , 
dER m x tua J v . dEn, 



where the local DM density po = 0.3 GeV/cm and f m in — J m 2 A i R ■ dR/dEn includes contributions from both \ 



and its conjugate \ for they have the same cross sections. 

In the non-relativistic limit, the differential cross section can be Maclaurin expanded in powers of v r . The two 
most important contributions are 



da 



-d{a-} + d{a+} 



with v r independent coefficients denoted by brackets. (For example, in Eq. \T\ d{a-} is the coefficient of v^ 2 , and 
d{a+} — 0.) Usually, the first term is the relevant one (as in the EDM, CFF, or SI cases). However, in certain 
cases like MDM, the second term may compete due to the 1/Er enhancement from the low energy virtual photon 
propagator. On integrating, we find 



dR po 1 
dEa m x vha 

where the dimensioless integrals are defined by 
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with the shorthand v u — min(t> m i n + de,ii csc ), Vd — min(t; m i n — ve,v C ec)- Note that /_ = for v m { n > v CEC + ve- 



3 Magnetic dipole moment of dark matter 

In the static limit, the DM magnetic moment can only couple to the nuclear magnetic moment via 
the induced magnetic field. The non-relativistic effective Hamiltonian of a magnetic moment of a 
particle with spin S subject to a magnetic field B is, in our convention 

H cS = -fiB-S/S, 

which agrees with the standard form for a spin i particle, i.e., —fiB • <x, where \x is dimensionful. 
Since B = V x A, we identify the curl Vx with the momentum transfer qx . The direct scattering 
off the nucleus of spin / via the interaction between the nuclear magnetic moments nz,A an d the 
magnetic moment fi x of the DM particle is described by 

The curl substitution has been applied to the DM vertex and the nuclear vertex which are linked 
by the virtual photon exchange factor 1/g 2 . The spin-averaged differental cross section is 

A I AJ\ 2 2 2 S + lI+ldq 2 22 

da M dm UN) = ~^x^z,A-^g 3i~4^2 ' M ^ 9 » ' &' 

The magnetic nuclear form factor Gm has been included. The above purely magnetic description 
ignores an important nuclear charge Z 2 effect which is suppressed by v 2 but enhanced by 1/Er. 
Both the magnetic and electric effects are comparable at direct search experiments with low recoil 
energy thresholds. 

To account for the Z 2 effect, we need to treat the convection current of the nucleus Z(p% A + p'% A ) 
(where p^ z A and p% A are the incoming and outgoing momenta, respectively) with a Dirac trace 
calculation. Then we include the above contribution from the nuclear magnetic moment associated 
with its spin /. The leading contributions for low recoil energy and low relative velocity give 
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The two nuclear form factors Ge and Gm are normalized so that G#(0) = 1 and Gm(0) = lj 
the charge Z and the moment fJ,z,A have been factored out. In our numerical analysis we make the 
simplifying assumption that the two nuclear form factors Ge and Gm are approximately equal, 
although they could be slightly different. 

Several studies only include the nuclear spin and ignore the effect of the nuclear charge, e.g., 
Ref. [6] . Others inappropriately scale the result for a Dirac point charge by Z 2 [TJ [8] . Our result 
reduces to that of Ref. [8] in the limit that a nucleus becomes a Dirac point charge, Z — > 1, 
Hz, A I '27I *" 1) rn A —> m p> and I — > 1/2, and has subsequently been reproduced in an erratum to 



Ref. [9]. It is worth pointing out that a proton with a significant anomalous magnetic moment is 
not faithfully represented by a Dirac point charge. Note that the A 2 factor in the result of Ref. |10j 
must be replaced by (mA/m p ) 2 to agree with our expression, and for quantitative accuracy. 

4 Electric charge form factor of dark matter 

A neutral non-self-conjugate dark matter particle can couple directly to the photon via its charge 
form factor (CFF) G x (q 2 ); charge neutrality only implies G x (0) = 0. At low momentum transfer, 
it is conveniently parameterized by a cutoff Aqff'- 

G x {q 2 ) = q 2 /A 2 CFF . 

The amplitude for DM-nucleus scattering via the form factor is given by 

M = G x (q 2 )(e 2 /q 2 )ZG E (q 2 ), 

which is clearly spin-independent, so that the differental cross section is 

e*\ZG E (q 2 )\ 2 dq 2 
d "CFF( X N)= ^^ ^- (5) 

Employing the cutoff Aqff translates into a useful cross section for DM-proton scattering: 

M 1 mlm 2 / e 2 \ 2 / 1 TeV \ 4 m\ , 2 

a%> = — - p x ^ — ~ *— _ x 0.92 x 10~ 42 cm 2 



T (P) = 1 m P m x ( <? V^ ATeV y m x 
CFF it (m x + m p ) 2 \A 2 CFF ) V a cff J (m P + m x ) 2 



which is relevant to the DM capture rate in the Sun. The proton form factor is essentially flat and 
ignored during the integration. 

In comparison, the SI scattering amplitude, 

Msi = 2[f p Z + f n (A - Z)]G E (q 2 ) , 

(where we have used the electric form factor to approximate the nuclear size effect), yields the 
differential DM-nucleus cross section, 

and the DM-proton SI cross section, 

Jp) _ ( Z Y Ap) _ 4 m 2 p m 2 x 
a si~[ A ) a CFF-- {rnx + mp) 2h- 

We therefore have the correspondence, 

e 2 

Z ^ 2[f p Z + f n (A - Z)\ . 
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If we assume f p = f n = A S7 for simplicity, we can further simply to the correspondence, 
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A S7 2 \ A J A CFF 

Since 2Z ~ A for most nuclei, the correspondence becomes 

2 A 

A SI O -Acff ~ 6.6 A C FF ■ 
e 

5 Analysis 

For the MDM of DM, we define the dimensionless gyromagnetic ratio or g- factor g x by 

&x s = ze: ■ 

2m x 

It is useful to compare d x and \x x in terms of the effective cutoff scales Aedm and Amdm defined 
by 



d x = 7 and ^x = Sx S -< 



e 



Aedm ' x 2m x A MDM 

Then, 

Aedm in _ 00 10" 21 e-cm Amdm 2 

19733 x and 



1 TeV d x 1 TeV g x S 1 TeV 

(Note that the cutoff scales are arbitrarily defined quantities, related to compositeness or short 
distance physics. They are defined in order to facilitate comparison and are not to be interpreted 
as new physics scales.) Direct DM search experiments are currently sensitive to d x ~ 10~ 21 e-cm 
which translates into a sensitivity to Aedm of almost 20 PeV. On the other hand, experiments are 
sensitive to g x S ~ 0.001 or Amdm of order 10 TeV for m x < 10 GeV. The difference arises because 
the MDM contribution to direct searches is more non-relativistically suppressed. Due to the higher 
dimensionality of the DM charge form factor, the corresponding scale Acff can be probed to about 
a few 100 GeV. 

We now consider the consequences of these electromagnetic properties of DM for direct searcheso 
Since the DM velocity cannot exceed v esc , which we take to be 600 km/s, a DM particle lighter than 
7.2 GeV cannot give a nuclear recoil energy above 10 keV, which was the threshold cut adopted 
by the CDMS II collaboration [12]. Although CDMS II has accumulated more kg-days of data 
than CoGeNT |13| , the lower threshold energy of CoGeNT permits sensitivity to masses as low as 
3 GeV. Solely for the purpose of illustration, we focus on the CoGeNT anomaly. 



4 The relevant nuclear parameters for the CoGeNT, CDMS and XENON experiments are |llj . 



73^^. r _ 9 



Ge: I = f , hz.a = -0.8795^-, / = 7.73% 



where / is the natural abundance of the isotope. 



IfXe: I = i, vz, A = -0.778^, / = 26.44% 
ifXe: / = §, M z,A = 0.692^, / = 21.18% , 
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Figure 1: A DM particle with m x = 7 GeV and an EDM, MDM or CFF reproduces CoGeNT data 
from a 56-day run with 0.33 kg of germanium. Note that the shape of the CFF curve is identical 
to that for spin-independent scattering. 



In Fig. 1, we show that a 7 GeV DM particle with an EDM of d x = 10 20 e-cm (equivalently 
Kedm = 1.97 PeV), or an MDM with g x S = 0.00454 (equivalently A M dm = 3.09 TeV), or a CFF 
with Acff = 187 GeV (corresponding to Oq FF = 5.8 x 10 -40 cm 2 or agj = 1.1 x 10~ 40 cm 2 ) easily 
reproduces the CoGeNT event excess below 2 keVee; we employed an energy-dependent quenching 
factor, E'(keVee) = 0.19935 E^keV) 1 ' 1204 , which converts the total nuclear recoil energy En to the 
energy detected by the experiment (in the form of ionization, scintillation or heat) with units of 
equivalent electron energy (keVee). The corresponding x 2 values for the 3-parameter fit to 30 data 
points between 0.45 keVee and 1.9 keVee are Xedm = ^6, Xmdm = ^2 and Xc FF = 20, where [2] 



X 2 = £ 2(iVf 

i with Nf XI V0 



N exp 

1 i with Nf xp =0 



2N 



tii 



and the Nf include a background contribution that is modeled as a linear combination of a constant 
term and a sum of two weighted Gaussian distributions that describe the peaks from the decay 
of 65 Zn and 68 Ge via L-shell electron capture |13j . Two of the fit parameters fix the background 
and the third parameter normalizes the signal. In all cases, the constant background is about 3 
events/bin which explains the the data between 1.5 keVee and 3.2 keVee satisfactorily. 

Since in all cases the DM-nucleus scattering is dominantly SI , the most stringent constraints in 
this low m x range are obtained from the XENON10 experiment due its low energy threshold |15| 
(and perhaps CDMS-II silicon detector data |16] which may have an underestimated energy calibra- 
tion uncertainty). We assert that the DM candidates of Fig. 1 are consistent with the XENON10 
data based on the following: (1) The shape of the CFF recoil energy distribution is identical to 
that for SI scattering, and depending on analysis details, SI scattering of a 7 GeV DM particle 
with (jgj ~ 10~ 40 cm 2 is either compatible with [T7] or marginally excluded by |18] XENON10 
data, (2) the distributions of Fig. 1 are almost identical, and (3) we checked that for m x = 7 GeV, 
the difference in the recoil energy distributions for CFF scattering on Ge and Xe is larger than for 
scattering via an EDM or MDM. Taking these facts together, we infer that the DM candidates of 
Fig. 1 are not in conflict with XENON10 data. 
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